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Abstract. A superintegrable finite model of the quantum isotropic oscillator in two 
dimensions is introduced. It is defined on a uniform lattice of triangular shape. The 
constants of the motion for the model form an SU (2) symmetry algebra. It is found that 
the dynamical difference eigenvalue equation can be written in terms of creation and 
annihilation operators. The wavefunctions of the Hamiltonian are expressed in terms of 
two known families of bivariate Krawtchouk polynomials; those of Rahman and those 
of Tratnik. These polynomials form bases for SU (2) irreducible representations. It is 
further shown that the pair of eigenvalue equations for each of these families are related 
to each other by an SU(2) automorphism. A finite model of the anisotropic oscillator 
that has wavefunctions expressed in terms of the same Rahman polynomials is also 
introduced. In the continuum limit, when the number of grid points goes to infinity, 
standard two-dimensional harmonic oscillators are obtained. The analysis provides the 
N — > oo limit of the bivariate Krawtchouk polynomials as a product of one-variable 
Hermite polynomials. 
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1. Introduction 

There is a considerable body of knowledge on superintegrable models and remarkably, 
the recent years have witnessed significant advances in their classification. The harmonic 
oscillator is the prototype of superintegrable models. These are Hamiltonian systems 
which admit a number of conserved quantities greater than the dimension d in which they 
are defined. They are called maximally superintegrable when this number is 2d — 1. In 
quantum mechanics, these constants of motion are typically required to be algebraically 
independent. To a large extent, the documented cases describe continuous systems and 
little is known in the discrete realm. It is the purpose of the present paper to present a 
finite analog of the quantum harmonic oscillator in two dimensions which is maximally 
superintegrable. 

There is by now a number of finite models of the quantum oscillator in one 
dimension. The most studied is based on the SU(2) algebra pp. Generalizations 
that model the parabosonic oscillator have been recently developed with extended 
algebras [131 HU Let us recall the main features of the system associated with 
the 577(2) algebra whose generators Jo, J± satisfy [Jo, J±] = ±J± an d [J+, J_] = 2J 
and whose standard N + 1-dimensional irreducible representations may be taken to 
have basis vectors \N, n), n — 0, . . . , N, defined by Jo|iV, n) — (n — N/2)\N, n) with N 
a non-negative integer related to the value of the Casimir operator. Upon making the 
identification 

H = Jo + j + l, Q = \{J + + J-), P= Z -(J + -J-), (1) 

the SU (2) commutation relations are interpreted as the Heisenberg equations of motion 
of the oscillator: 

[H,Q] = -iP, [H,P)=iQ. (2) 

The wavefunctions of this finite oscillator are, naturally, the overlaps between the 
eigenstates of the position operator Q and those of H, ie. {\N, n)}; they are hence 
given by rotation matrix elements known to be expressible in terms of the Krawtchouk 
polynomials [IE] K n (x, |, TV) with x = N/2 + q and q G {—N/2, N/2}, the eigenvalues 
of Q. The Krawtchouk polynomials are given in terms of a hypergeometric series 



N 



K n (x;p;N) = (_A0„ g !Z!iMZ±2i (j. j (3) 

and are orthogonal with respect to the binomial distribution 

N 



x 



p*(l-p)"- x , x = 0,-.-,N. (4) 



This entails a realization of the observables of this finite oscillator in terms of finite 
difference operators. It is known that S77(2) contracts to the Heisenberg algebra. Hence, 
when N — > oo and with proper rescaling, the description of the standard (continuous) 
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harmonic oscillator is recovered [T] as the Krawtchouk polynomials tend to the Hermite 
polynomials H n (x) in the limit iV — > oo. 

The Krawtchouk (or SU(2)) one-dimensional finite oscillator has been used to build 
two-dimensional models. A first approach was to take the direct product of two finite 
one-dimensional oscillators to obtain a two-dimensional system defined on a square grid 
of points and with SU(2) x 577(2) as dynamical algebra. A second one [I] was to use the 
isomorphism between SO (4) and 577(2) ®SU(2) to describe this finite two-dimensional 
oscillator in terms of discrete radial and angular coordinates. As remarked in Refs. 
[H IH E] , these models do not possess the SU (2) symmetry algebra of the continuous 
two-dimensional quantum oscillator. This symmetry algebra (not to be confused with 
the dynamical algebras) is generated by the three constants of motion which make the 
continuous harmonic oscillator maximally superintegrable. The two-dimensional finite 
models constructed so far do not therefore share that property. We here introduce a 
system that has such an SU(2) symmetry. The model will be defined on a grid of 
triangular shape. 

Finite oscillator models have various uses e.g. in optical image processing and in 
signal analysis where only a finite number of eigenvalues exist [2j [3]. Finite planar 
oscillators can be employed in particular to describe waveguides and pixellated screens. 
It is expected that preserving, at the discrete level, all the symmetries of the continuous 
oscillator model would prove advantageous in these applications. 

It is natural to think that two-variable generalizations of the Krawtchouk 
polynomials could provide an underpinning for interesting two-dimensional oscillator 
models. As a matter of fact, two families of such bivariate polynomials have been 
identified. Some 20 years ago, Tratnik introduced multivariable Racah polynomials 
depending on two parameters, that encompass multivariable Krawtchouk polynomials in 
their simplest case. Geronimo and Iliev j6] subsequently showed that these polynomials 
are bispectral, i.e. that they obey, in the bivariate case, a pair of difference equations 
in addition to the recurrence relations associated with their orthogonality. Griffiths 
[7] in the early 70's and more recently Hoare and Rahman [11] discussed a second 
family of bivariate Krawtchouk polynomials with four parameters in connection with 
a probabilistic model. (See also Ref. [18J.) These were later extended to an arbitrary 
number of variables [9] and shown to be also bispectral [10J. In contradistinction with 
the Tratnik version, the Krawtchouk polynomials of Rahman satisfy a nearest neighbor 
five-term recurrence relation [8]. They have been used to construct a spin lattice with 
remarkable quantum state transfer properties [T7| . 

In view of the bispectrality of the two-variable Krawtchouk polynomials, it is 
of interest to inquire what is the quantum mechanics that their difference equations 
describe and to find the corresponding continuous dynamics in the limit N — > oo. 
This last question involves obtaining the iV — > oo limit of the bivariate Krawtchouk 
polynomials, an issue that seems not to have been resolved so far. 

The Rahman polynomials lead to the most direct interpretation. Remarkably, 
their iV — > oo limit is found to be the product of two Hermite polynomials. This 
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is obtained by observing that their difference equations factorize in terms of a pair 
of simple creation and annihilation operators which obey Heisenberg commutation 
relations. A finite model for the two-dimensional isotropic (or anisotropic) oscillator 
can then be constructed in the standard fashion. In terms of these creation and 
annihilation operators, the constants of motion take the same form as in the continuous 
situation and the symmetry algebra is therefore preserved in the discretization process. 
For the isotropic case, the Rahman polynomials are shown to form a basis for the 
irreducible representation spaces of SU(2) They correspond in the continuum limit to 
wave functions that are separated in Cartesian coordinates. In view of the eigenvalue 
equations that they obey [6J, the Tratnik polynomials are seen to correspond to 
separation of variables in coordinates that are rotated relative to the ones associated 
with the Rahman polynomials. This is why the Tratnik polynomials depend essentially 
on one less parameter than those of Rahman; the additional parameter has been removed 
by the rotation. It should be stressed that the oscillator models based on these bivariate 
Krawtchouk polynomials are defined on a uniform grid of triangular shape. 

The remainder of the paper is organized as follows. In section 2, we present an 
SU (2) invariant finite isotropic oscillator model, whose eigenfunctions are given in terms 
of the two- variable Krawtchouk polynomials of Hoare and Rahman, and we describe its 
continuum limit. Also, a finite anisotropic oscillator is introduced which is seen to 
be superintegrable for rational frequencies of the coupling constants. In section 3, we 
show how the two-variable Krawtchouk polynomials of Tratnik can also be realized as 
eigenfunctions of the same isotropic Hamiltonian, though associated with a different 
integral of motion. Conclusions and an Appendix on the large N limit of the trinomial 
distribution follow. 



2. Finite oscillator model and the two-variable Krawtchouk polynomials of 
Hoare and Rahman 

We first begin with the two-variable Krawtchouk polynomials introduced by Hoare and 
Rahman [TT] which are a special case of the Aomoto-Gelfand hypergeometric function 

m 

m ' n ~ ^ i\j\kW.(-N) i+j+k+e W 

o<i+j+k+e<N J v ji+j+K+z 

with parameters 

u _ (Pl+P2)(Pl+Pa) _ (Pl+P2){P2+Pi) 

1 piEp. ' 2 P2EPi (p.) 

y _ Ql+P3)(P3+P4) _ (P2+P4)(P3+P4) V U / 

It should be noted that the polynomials depend essentially on three independent 
parameters. This fact can be observed in terms of the Pi by writing 

u = (1+P2/Pl)(l+P3/Pl) / 7 v 

1 1 + P2/P1 +P3/P1 +P4/P1 

and equivalently for u 2 ,Vi and t> 2 . 
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With m and n nonnegative integers such that < m+n < N, the polynomials K^ n 
are defined on a uniform lattice of a triangular shape, i.e. {(x, y) G Z> | < x+y < N} 
and are orthogonal with respect to the two- variable generalization of the weight (jlj) given 
by the trinomial distribution 

u(x, y)=( N ) «(1 -m- V2) N ~ x - y , (8) 



with 



x,y 

(Pl+P2)(Pl + Ps)(P2+Pa)' 



= PsPa(Pi +P2 +P3 + Pa) , g x 

2 (P2 +P4>{P3 + Pa){Pi + P3)' 
It should be noted that the constants Ui,Vi,t]i satisfy the following three functional 
relations [8]: 

U1V1 + vim = !> 
U2V1 + V2V2 = 1, 

u 1 u 2 T] 1 + v 1 v 2 r]2 = l, (10) 

which also imply that the polynomials K^ n depend essentially on three parameters. 

One of the specific features which the polynomials K^ n possess is that they obey 
a nearest neighbor difference equation in x, y [10J (and also some contiguity relations in 
m, n because of their duality): 

[(pi + p 3 )m - (p 2 + p 4 )n]K^ n 

— (N — X — y) J^»(P8+P4)Eft^ _ PsP^+P2)J2Pi A \ K N 
{ P1+P2 P3+P4 J 

Pl + P2 P3 + P4 

where A is a difference operator defined by A±/ 1): y a) .../(Zi, Z2, • ■ •) = f(h ± 1, ^2 ± 1, • ■ •) — 
f(h,h,"')- Recently, it was shown [121 HO] that the polynomials K^ n additionally 
satisfy a seven point difference equation, which is equivalent to the fact that the 
polynomials K^ n are eigenfunctions of the following difference operators: 

a n = ( N - x -y)pmT,Pi f P2 A _ Pi A 

1 (Pi +P3){PlPi -P2P3) \Pl +P2 X P3+P4 V 

P1P4-P2P3 ( P3 A Pi A 

+ 7 w r x A_ x — y A_ 

(Pi + P2) (P3 + Pa) \ P1+P3 Pi+ P3 



PsPA a PlP2 
' /\ x — y 

(Pl + Ps) (?3 + P4) (Pi + P2) (Pi + P3) 



a jv = (N-x- y)p 2 Pi E P» / Pi A ] P3 A 

" " (P2 +P 4 )(PlP4 -P2P3) V P1+P2 1 P3+P4 1 
P1P4-P2P3 / P4 A . P2 A 



(Pl + P2) (P3 + Pa) V P2 + P4 P2 + P4 

- x P?P* A _ y ^ A (12) 

(p 2 + Pa) (p 3 + p 4 ) (Pi + P2) (P2 + Pa) *' v ' V 7 
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The eigenvalue equations read: 



lv l Jx m,n 
iv 2 -"-m.n 



mK N 



nK, 



in , 1 1 ' 



(13) 

As it will later become important, the sum A^ + A^ depends only on the two parameters 
r/i and 772- 

A different combination of the operators flT2l) eliminates the non-nearest neighbor 
terms involving A x _ yi ^~x,y and the difference equation ( TTTT) recovered. Moreover, one 
can easily find that p 2 (pi + Ps)^i + Pi(P2 + Pa)^2 an d PiiPi + P3)Af + Ps(p2 + Pa)^2 
give the "recurrence" relations for the polynomials K^ m : 

(pa(pi + Ps)m + pi(p2 + p 4 )n)K^ n (x, y) 

{PiPa -P2P3) 2 



(N 



x 



v) A- 



+ x 



Pi +P2 
P3Pii.Pl +P2) 



X 



(Pi +P2)(P3+P4, 



-A 



P3 + Pi 

(p 4 (pi + P3)m + p 3 (p 2 + P4)n)K^ n (x, y) 



A- X ,y + yPlP2^x-y 



K 



N 



(N-x-y) 



pm E p* 

P3 +P4 



(PlP4 - P2P3)' 



(Pi +P2)(P3 + Pi 



■A. 



x ,—y 



(14) 
1 at 



Pi + P2 

If we set i^ m>n = when x or y is a negative integer and given that K n 
x = y = from (j5J), we can determine all the K^ n (x,y) from (TT4"j) . In that sense, the 
operators A^ and A^ fully characterize the Rahman polynomials K^ n defined on the 
grid points such that 0<x + y<Na.s their joint eigenfunctions. 



2.1. The SU(2) invariant isotropic Hamiltonian 

In this section, we first consider the operators defined by 

AR) _ PlP2P3P4T,Pi 



A 

A (L) 



(Pi + Pz)(PlPi - P2P3) \P4(Pl + P2 

P1P2P3PiY.Pt ( - 1 



(P2 + Pi){PlPi ~ P2P3) \P3(Pl + P2 
Hi.. A ! PlP4 ~ P2P3 

' '' P3 " 

T7 1 + — T- 



1 

-A x . -. rA„ 

1 P2(P3+P4) 

"A, + — — -Ay 

1 Pl{P3+P4) 



A, 

A (L,N) 



P1+P2+P3+ Pi \Pl 

pm - P2P3 ( -x r ^ 

P1+P2+P3+P4 V P2 



P4 



+ (N + 1 - x - y) 
+ (N + l-x-y), (15) 



where T is the shift operator defined by Tif(l) = f(l + 1). One can easily verify that 
these operators satisfy a shifted form of the Heisenberg algebra relations: 



A (i) A (i,N) _ A (i,N-i) A (i) 



R,L 



A iR) A (L) 



A [L) A {R) 
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A (R,N+1) A (L,N) _ A (L,N+1) A (R,N) _ Q 
A (R) A (L,N) _ A (L,N-1) A (R) = Q 

A (W) A (L) _ + A (L) A (R,N) = Q ' (16) 

and provide the factorizations of the operators and : 

Af = Ap N -Q A m ) A % = A f N - 1] A^. (17) 
Let us consider the polynomials A_ K^ n . It is straightforward to check from ffT6j) that 

Af- 1 (AW< B ) = (m-l)^< B , 

A?- 1 ^^)^^^. (18) 

This means that A ( * ] K^ n are also the Rahman polynomials K^_\ n up to a 
multiplicative constant, which is determined by evaluating A^K^ n at (x, y) = (0,0). 
We can thus obtain the following four relations: 

A( R )fr N =—K N ~ 1 A {L) K N =—K N - 1 

— m,n pj in—l,ni — m,n j\j m.n— 1' 

Af N) Kl n = (N + A + ,N) Kn = (N + l)K»f +1 . (19) 

Hence, the operators (1151) provide ladder operators for the Rahman polynomials (J5J). 
Now let us go back to the relations (fT6|) . From these, it is clear that the Hamiltonian 

f) JSO = Af + A^ (20) 

admits the following integrals of the motion 



Jx = - (Ai R ' N ^Ai L) + A^ N ~ l) A w 



1 

' X = 2 

J Y = l - (Af N - x) A^ - A$f-*A™ 

1 f ,(R,N-1) A (R) AL,N-X) AL) 



j z = l i^A^-^A^' - A^' iy - l) A^j . (21) 
The conserved quantities ( I2T1) form a basis for an SU(2) algebra as they satisfy 

[Jx,Jy]=Uz, [Jy,Jz]=iJx, [Jz,Jx]=iJy (22) 

and hence the Hamiltonian \) iso of ( 120]) is SU(2) invariant. It is readily seen that the 
Casimir operator 

Q = J\ + 4 + Jf (23) 
takes the form 

Q=j(j + l), J = ^- (24) 

It hence follows that the Rahman polynomials fl5]) with fixed j = (m + n)/2 form a basis 
for the 2j + 1-dimensional irreducible representations of SU(2). They diagonalize the 
operators A^, A^ or equivalently the commuting pair t) iso and Jz with eigenvalues 

t, iso K» >n = (m + n)K» n , J z K^ n = \{m- n)K^ n . (25) 
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2.2. An anisotropic Hamiltonian 

Anisotropic oscillators can also be constructed in the standard way using the ladder 
operators ( TT5l) and the "number" operators and factorized in ( TT71) . Systems 
with 

f) amso = u; 2 Af (26) 

as Hamiltonians will be integrable for arbitrary u\ and u 2 and superintegrable when the 
ratio of the frequencies is rational. In each case, the eigenf unctions for the pair \) an iso 
and either A^ or A^ are the Rahman polynomials. It is important to note that whereas 
f)i SO depends essentially only on two parameters r\\ and 772, the anisotropic Hamiltonian 
depends on all three functionally independent parameters. 

An interesting observation in this context is that the spin lattice Hamiltonian 
introduced in [17] is generically associated with an anisotropic oscillator. In fact, a 
Hamiltonian of the form ( 126]) will have nearest-neighbor interactions only if the relation 

= (27) 

\U 2 J P2+P4 

is verified. 



2.3. Continuum limit 

In this subsection, we consider the N — > 00 limit of the two-variable Krawtchouk 
polynomials (JSJ). The limit will be obtained by the following change of variables 

x = Nt]i + VN (cis + c 2 t) , y = Nrj 2 + Vn (c 3 s + c A t) , (28) 

with 



Cl 



c 2 



C3 



-P2 /2(pi +p 2 + p 3 +Vi)vm 



(Pi + P2) (Pi + Pa) V P2+P4 



Pi /2(pi +P2 +P3 +P4)P2PA 



(Pi + P2) (P2 + P4) V Pi + P3 



P4 /2(pi +p 2 +P3 +P4)P2P4 



(P3 +P4)(Pl +P3) V P1+P3 



^ = -P3 / 2(P1 +P2 +P3 +P 4 )P2P4 ^ 

4 (P3 +P4)(P2 +P4) V P1+P3 

Note that as A" — >■ 00, the range of s and t becomes the whole real line. 
The operators A^, A^ have the following limits 

lim Af = -\d 2 s + sd s , lim A^ = -\d 2 t + td t . (30) 

N— too 2 iV— >oo 2 

Thus, in the limit A" — > 00, the Hamiltonian i)i SO given in (120]) tends to the Hamiltonian 
of a two-dimensional oscillator 

lim f) iso = -i («9 S 2 + d 2 t ) + s<9 s + td t (31) 

TV— >oo 2 
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which has been conjugated by the ground state and with the coupling constants absorbed 
into the variables. 

While the operators and do not change the quantum numbers m, n and N, 
the ladder operators do and so the the normalization of the polynomials (see below) 
affects the limit of these operators. In order to find the proper choice of normalization, 
note that in the large A" limit the trinomial distribution tends to the two variable 



Gaussian (for a proof see Appendix A) with the following coefficients 



u(x,y) 



N 



z,V J 271-^^771772(1-772-772) 
Note that the constants in the change of variables (j29j) satisfy 
C1C4 - c 2 c 3 = 2a/t7i77 2 (1 - 771 - 772) 



+ N 



.(32) 



and so the unit of and y are shifted by 1, is (2^^/771772(1 — 771 — 772)) 1 . Thus, 

1 



OO POO 



00 J — OO 



lim F(x(s,y),y(s,t)) 

,N->oc 



lim V F(x,y)uj(x,y) = - 

/V-S'OO Z ' TT 

0<x+y<N 

The normalization of the Rahman polynomials given in j8] is 

0<x+y<N 

with 

mi\ni\(N — mi — ni)\(pip A — p 2 p 3 ) 2( - mi+ni,) 



1 dsdt. (33) 



(34) 



jni,n 2 
mi ,rri2 



"mi,m.2"ni, 



,na An(pip 3 (p 2 + P4)) mi (p 2 p 4 (Pi +p 3 )) ni (EPA ; ) mi+ " 1 

This suggests the following redefinition: 



.(35) 



K N ( x v) = a N K N (x v) 



iV 



Q', 



V 



AH {pip 3 (p2 + Pa))™ {p2Pi{pi + Pz)) n 



(piPi — P2Pz) m+n V (A^ — 771 — 77)! 

With this normalization, the orthogonality relation becomes 

£ y)K% >n K% >n , = 2 m+n n\m\5 nin ,5 m , m/ , 



0<x+y<N 

which implies, in the limit 

1 /"OO />oo 
7T 



00 >/ — OO 



lim K, 



N 



7V-5>oo 



dsdt = 2 m+n n\m\. 



(36) 



(37) 



(38) 



Thus, the functions limjv-j.00 K^ n satisfy the eigenvalue equation 
1 



§3 + sd. 



(39) 



lim K% ) = m ( lim AT^ 

viV— >-oo ' J \N^oo ' / 

and a similar one in t. They are moreover square integrable in the plane with respect 
to the measure e~ s2 ~ t2 and are hence polynomial. The normalization constant has been 
chosen to agree with a product of Hermite polynomials and so the limit satisfies 

lim K» n = eH m (s)H n (t), e 2 = 1. (40) 

N->oo ' 
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Using the proper normalization, the action of the ladder operators becomes 



yipmiPi + P2 + P3 + Pa)(P2 + Pa) . (m,N) f>N _ f>N+l 
[PiPa — PiPz)v N + 1 



^4— K mn — TnK m _i n} (41) 

y / 2p 1 p 3 (p 1 + p 2 + p 3 + Pa)(P2 + Pa) 

and similarly for the operators A\. As expected, the limit of these scaled operators 

becomes the ladder operators for the Hermite polynomials 

s Um _ V2PMP1+P2 +PS+ Pa) fa + Pa) Air , N) = 2s _ 



a m = lim -^N(PiPa-P2Pz) ^n) = \ Q 



2piP 3 (Pi + P2 + P3 + Pa) (P2 + Pa) 2 
\I^P2Pa{p\ + P2 + Ps + Pa) (pi + Pa) ,(n,iV) 



a" = lim _ v-™^_^ ; ^ , ^w ^^.; = 2i _ ^ 

a" = lim -^(P^-PzPs) A jn,m = 1 ^ 

n^oo y/2p 2 pA(pi +P2+P3+ Pa)(Pi + Pa) 2 
Since these ladder operators correspond exactly to those of Hermite polynomials, the 
sign in f HOj) does not depend on either m or n. Furthermore, since K^ = H (s)H (t) = 
1, the sign in (jlDj) is exactly e — 1. 



3. Two-variable Krawtchouk polynomials of Tratnik 

A different set of commuting operators in the enveloping algebra of ( fl5j) lead to the 
version of two- variable Krawtchouk polynomials defined by Tratnik [19, 20J. Beginning 
with the observation that the isotropic Hamiltonian f)j so depends only on the two 
parameters % and rj 2 , it is also possible to construct another operator written in terms 
of the creation and annihilation operators ([151) . given ffTTl) . that depends only on these 
two parameters. The set of commuting operators is given by 



Co 



p 2 p,A^ N) A^ + p m Af N) A^ + PlP4 Af + P2P3 A? 



where 



(Pa + Pa) (Pi + Pz) (P2 + Pa) (Pi + Pz) 



PiPa P2P3 



,(43) 



,P2+P4 P1+P3, 

Thus, the joint eigenfunctions of these operators will depend essentially only on two 
parameters. These eigenfunctions are the two-variable Krawtchouk polynomials of 
Tratnik as the operators L\ and C 2 are exactly those given in Geronimo and Illiev [B]. 
Since C 2 is a difference operator in y alone, the two-variable Krawtchouk polynomials 
can be written as imbricated ordinary Krawtchouk polynomials 

K 2 (n; x, y; p ; N) = ^(x,Pi,N-n 2 ) ^ / . fa ^ \ {u) 

{— iV )n 1 +n 2 V 1 — Pi / 
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where 

= i = pmYsPi 

1 1 (Pl +P2)(Pl +P3)(P2 +P4)' 

p 2 = m = pmY^Pi / 45 n 

2 2 (P3 +P4)(Pl +P3)(P2 + Pi)' 

The eigenvalue equations are 

dK 2 (n; x, y; p; AT) = (m + n 2 )K 2 (n; x, y; p; iV), (46) 
£ 2 iT 2 (n; x, y; p; AT) = n 2 K 2 (n; x, y; p; A/ - ). (47) 

While both the Tratnik and Rahman polynomials diagonalize the isotropic Hamiltonian, 
the Rahman polynomials depend on an additional parameter. The reason for this is that 
the two families of polynomials each diagonalize a pair of operators which are related 
by an SU(2) automorphism. 

Recall from Section |2] that the Rahman polynomials are eigenf unctions for the 
isotropic oscillator Hamiltonian f) iso and the operator Jz- The Tratnik polynomials are 
also eigenfunctions for t) iso as well as the operator 

t=±C 1 -£ 2 , (48) 
with eigenvalues 

lK 2 [n- x, y; p; N) = i(m - n 2 )K 2 {n; x, y; p; N). (49) 

Because of their eigenvalues, it is natural to expect that the operators Jz and t are 
related and indeed t can be expressed as 

t = aJ x + bJ y + cJ Z} a 2 + b 2 + c 2 = 1 (50) 



with 



P2 +P4 Pl+ P3 
P2P4 PlP?, 



P2P. + ^P^ ] (51) 



P2 +P4 Pl+ PZ 



(52) 



c = 2£ 2 -^- - 1. (53) 

Pl +P3 

This of course corresponds to the observation that t is obtained from Jz by a rotation. 

The Tratnik polynomials (jUj) also tend to products of Hermite polynomials in the 
continuum limit as their imbricated expression suggests. The limit of the eigenvalue 
operators for these polynomials takes the form 

lim C 1 = -l(d 2 + d 2 )+sd s + td t , 

iV-KX) Z 



lim C 2 

N^oo 



PiPa d 2 p 2 p 3 d 2 



2(Pi + P2)(P2 + Pa) ds 2 2(pi + p 2 )(pi + Ps) dt 2 

Pms I P1P2P3P4 A d 



(Pl + P2){.P2 + Pi) V (Pl + P2) 2 (Pl + P3)(P2 + Pi) J ds 
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+ 




(Pi +P2)(Pl +P3) 



(pi +V2) 2 {pi+ P3)(P2 + Pi) dtds 



P2P3t 



P1P2P3P4 



+ S 



(Pi +P2) 2 (Pl +P3)(P2 +P4) 
d 2 1 



PlP2P-ZPi 



) 



dt 



(54) 



As expected, the operator £2 can be rotated into a linear combination of two eigenvalue 
operators for Hermite polynomials in the new coordinates. 

4. Concluding Remarks 

We have introduced and described a finite oscillator model in two dimensions which 
has the same symmetry as its continuum limit, the standard two-dimensional isotropic 
harmonic oscillator. This finite model is constructed using eigenvalue operators for two- 
variable Krawtchouk polynomials. These operators were shown to factorize into ladder 
operators that obey Heisenberg-like commutation relations. The SU(2) symmetry 
algebra of the Hamiltonian is then obtained a la Schwinger using these raising and 
lowering operators and thus making the superintegrablity of the model manifest. As 
a consequence, the bivariate Krawtchouk polynomials were found to form bases for 
irreducible representations of SU(2). 

The Hamiltonian can be diagonalized by both the two-variable Krawtchouk 
polynomials of Rahman as well as those of Tratnik. In the two cases, the large N 
limits of the polynomials were obtained and seen to be products of Hermite polynomials. 
The two sets of polynomials are eigenfunctions of the Hamiltonian and of an integral 
of motion whose limit is a second-order integral associated with separation of variable 
in Cartesian coordinates. Additionally, we have shown how these additional integrals 
are related by an SU(2) automorphism thus providing a link between the two families 
of polynomials. We plan to return to this question in a future publication to provide 
explicit relations between Rahman and Tratnik polynomials. We have also provided 
models of anisotropic oscillators whose integrals can be constructed in the standard 
fashion from the ladder operators and whose wavefunctions are given by the Rahman 
polynomials without superfluous parameters. While we focused here, for simplicity, on 
two-dimensional problems, our study has a natural extension in arbitrary dimension; 
this shall also be the object of further investigation. 

Appendix A. Limit of the trinomial distribution 

Consider the limit of the trinomial distribution in the large N limit, with x and y 
renormalized as in (1281) . The proof follows the classical de Moivre-Laplace theorem for 
binomial distributions. For notational simplicity we denote 



x = Nrji + VNfi, 



y = Nr] 2 + VNu, 



(A.l) 



with 



fj, = C\S + c 2 i, v = c 3 s + c 4 t, 
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C1C4 - c 2 c 3 = 2^771772(1 - 771 - 772)- 

From Sterling's formula we obtain 

N 



13 
(A.2) 



W 



x,y 



N 



2n y xy(N — x — y) \Nrji 
AT — s — 7/ 



i + e>( 



Nrj 2 
1 , 



Making the change of variables (lA.ljl . allows for the approximation 

1 



A' 



xy(N-x-y) N^n^l - 7^ - r/ 2 ) 



(A.3) 



(A.4) 



and hence 



2^-^/7717/2(1 - r/i - 772) 



1 + 



1 + 



A^T/2 

1 



1 - 



'Nrn 
\l + v 



N(l - 7/1 - 7/2) 



2^1/7/17/2(1 - 7/i - 7/2) 



exp 



-x In ( 1 + 



A^T/i 



7/ In 1 + 



(JV - x - y) In ( 1 



+ 0{N~ 



y/Nil-rn-m)/. 



+ 0(N—2) 

Finally, a Taylor series expansion of the logarithmic terms gives 
N 



(A.5) 



7/i 7/|(l - 7/i - 7/2 



X,7/ 



,JV-x- 



9(1 — 770 — 779) ^ ^ 



2^^/7/17/2(1 - 7/2 - 7/2) 



(A.6) 
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